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Abstract 

r\ (%\ I 

We present measurements of the decays B —> D s D*~ , using data recorded by the BABAR 

detector in 1999 and 2000, consisting of 20.8 fb~ . The analysis is conducted with a partial re- 

(*)+ — 
construction technique, in which only the Dg and the soft pion from the D* decay are re- 
constructed. Prom the observed rates, we measure the branching fractions B(B° — > DfD*~) = 
(1.03 ± 0.14 ± 0.13 ± 0.26)% and B(B° D*+D*~) = (1.97 ± 0.15 ± 0.30 ± 0.49)%, where the 
first error is statistical, the second is systematic, and the third is the error due to the Df — > (J)tt + 
branching fraction uncertainty. From the B° —> D* + D*~ angular distributions, we measure the 
fraction of longitudinal polarization Tl/Y = (51.9 ± 5.0 ± 2.8)%, which is consistent with the 
theoretical prediction, based on factorization. These results are preliminary. 
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1 Introduction 



Precise knowledge of the branching fractions of exclusive B decay modes provides a test of the fac- 
torization approach used for the calculation of these branching fractions. Further tests are provided 
by measuring the polarization in decays of B mesons to vector-vector final states. Within current 
experimental sensitivities, these measurements are consistent with the factorization predictions for 
the final states D*p §, D*p(U50) @, and D*D* @. 

In this paper we present measurements of the branching fractions]] B(B° — > Ds + D*~). We 
also report the measurement of the D* + polarization in the decay B° —* D* + D*~, obtained from 
an angular analysis. These results provide increased precision tests of factorization. 

2 The BABAR Detector and Dataset 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage ring. 
An integrated luminosity of 20.8 fb _1 was recorded in 1999 and 2000 at the T(4S) resonance, 
corresponding to about 22.7 million produced BB pairs. 

Since a detailed description of the BABAR detector is presented in Ref . [||] , only the components 
of the detector most crucial to this analysis are briefly summarized below. Charged particles 
are reconstructed with a five-layer, double-sided silicon vertex tracker (SVT) and a 40-layer drift 
chamber (DCH) with a helium-based gas mixture, placed in a 1.5 T solenoidal field produced by a 
superconducting magnet. The charged particle momentum resolution is approximately (5pt/pt) 2 = 
(0.0013 pt) 2 + (0.0045) 2 , where px is given in GeV/c. The SVT, with a typical single-hit resolution 
of 10 pm, provides measurement of the impact parameters of charged particle tracks in both the 
plane transverse to the beam direction and along the beam. Charged particle types are identified 
from the ionization energy loss (dE/dx) measured in the DCH and SVT, and the Cherenkov 
radiation detected in a ring imaging Cherenkov device (DIRC). Photons are identified by a CsI(Tl) 
electromagnetic calorimeter (EMC) with an energy resolution a(E)/E = 0.023 • (E/GeV)^ 1 ^ © 
0.019. 

3 Method of Partial Reconstruction 

In reconstructing the decays B° — > Dg D*~ , with D*~ — > D°tt~, no attempt is made to recon- 
struct the D decays. Only the Dg and the soft 7r~ from the D*~ decay are detected. In this way, 

the candidate selection efficiency is higher by almost an order of magnitude than when performing 

(*)+ - 

full reconstruction of the final state. Given the four-momenta of the D s and tt and assuming 

that their origin is a B° — » Ds~*~ D*~ decay, the four-momentum of the B° may be calculated up 

(*)+ 

to an azimuthal angle eft about the Dl flight direction. This calculation also makes use of the 
total beam energy in the center-of-mass (CM) system and the masses of the B° and D*~ . Energy 
and momentum conservation then allows the calculation of the four- momentum of the D°, whose 
square yields the (^-dependent missing mass 

M miss = J {E hcam - E D (*)+ - E n ) 2 - (p B - P d m+ - p^-) 2 . (1) 

1 Reference to a specific decay channel or state also implies the charge conjugate decay or state. The notation 
Z)i*' + refers to either Df or Dl + . 
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In this analysis the missing mass is defined using an arbitrary choice for the angle (ft, such that the 
B° momentum ps makes the smallest possible angle with p^ and p n (*)+ in the CM frame. 



4 Event Selection 

For each event, we calculate the ratio of the second to the zeroth order Fox- Wolfram moment 0, 
using all charged tracks and neutral clusters in the event. This ratio is required to be less than 
0.35 in order to suppress continuum e + e~ — > qq events, where q = u,d, s, c. 



4.1 D+ and D*+ Candidate Selection 

S S 

We reconstruct Df mesons in the decay modes D+ -> <pir + , D+ -> K*°K+ and D+ -» K° S K + , with 
subsequent decays (ft —* K + K~ , K*° — ► K~ir + and K® — > ir + ir~ . These modes were selected since 
they offer the best combination of branching fraction, detection efficiency and signal-to-background 
ratio. The charged tracks are required to originate from within ±10 cm along the beam direction 
and ±1.5 cm in the transverse plane, and leave at least 12 hits in the DCH. 

Kaons are identified using dE/dx information from the SVT and DCH, and the Cherenkov 
angle and the number of photons measured with the DIRC. For each detector component d = 
{SVT, DCH, DIRC}, a likelihood L d (L^) is calculated given the kaon (pion) mass hypothesis. 
A charged particle is classified as a "loose" kaon if it satisfies jUl > 1 for at least one of the 
detector components. A "tight" kaon classification is made if the condition Yl d L^ /L^ > 1 is 
satisfied. 

Three charged tracks consistent with originating from a common vertex are combined to form 
a Df candidate. In the case of the decay Df — > (ftir + , two oppositely charged tracks must be 
identified as kaons with the loose criterion, with at least one of them also satisfying the tight 
criterion. No identification criteria are applied to the pion from the Df decay. The reconstructed 
invariant mass of the K + K~~ candidates must be within 8MeV/c 2 of the nominal (ft mass [||. In 
the decay — > 07r + , the (ft meson is polarized longitudinally, resulting in the kaons having a 
cos 2 9h distribution, where Oh is the angle between the K + and in the (ft rest frame. We require 
| cos 0h\ > 0.3, which retains 97% of the signal while rejecting about 30% of the background. With 
these requirements, the signal decay Df — > (ftir + and the Cabibbo-suppressed decay D + — > (ft-K + 
are readily observed (Fig. |l]a). 

In the reconstruction of the — > K*°K + mode, the K~tt + invariant mass is required to be 
within 65 MeV/ c 2 of the central K*° mass || . This wider window leads to a fraction of combinatorial 
background much larger than in the — > (ftn + mode. To reduce this background, we require 
|cos#ij| > 0.5. In addition, substantial background arises from the decays D + — > K*°tt + and 
D + — > K°tt + , which tends to peak around the nominal Df mass. This background is suppressed 
by requiring that the kaon daughter of the K*° satisfy the loose kaon identification criterion, and 
that the other kaon satisfy the tight criterion. Fig. |l|b shows the reconstructed K*°K + invariant 
mass. 

For the decay mode D+ -> K®K+, K° s — > 7T + 7r , the 7r + 7r invariant mass must be within 
15 MeV/ c 2 of the nominal K® mass, and the charged kaon is identified using the tight criterion. To 
improve the purity of the K® sample, we determine the angle a between the momentum and 
the flight direction defined by its decay vertex and the primary vertex of the event. We require 
cos a > 0.98 to reject the combinatorial background. The K®K + invariant mass distribution is 
shown in Fig. |l]c. 
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Figure 1: The invariant mass spectra for (a) 4"ir + , (b) K*°K + , (c) K®K + . The peaks are due to 
the D + (left) and Df (right) decays in this mode. The fit function is a single Gaussian for each 
peak, with the widths constrained to be equal, plus an exponential background. 
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Figure 2: Distribution of the mass difference AM = M D + — M D +. All the Df decay modes 
have been combined in this plot. The fit function is a Crystal Ball function for the signal plus a 
threshold function, as described in the text. 



The invariant mass M£> a of all Df candidate is required to be within three standard deviations 
( (J M Ds ) of the signal distribution peak M^ ak seen in the data. The standard deviations are a = 
5.23MeV/c 2 for </>7r+, 5.97MeV/c 2 for K*°K + , and 6.46MeV/c 2 for K°K+. 

All Df candidates satisfying all the above selection criteria are combined with photon candidates 
to form D* + — > Dfj candidates. The candidate photons are required to satisfy E y > 50MeV, 
where E^ is the photon energy in the laboratory frame, and E* > 110 MeV, where E* is the 
photon energy in the CM frame. When the photon candidate is combined with any other photon 
candidate in the event, the pair must not form a good ir° candidate, defined by a total CM energy 
£* 7 > 200 MeV and an invariant mass 115 < M 77 < 155 MeV/ c 2 . 

The distribution of the mass difference AM = M(Df"f) — M{Df) of events satisfying these 
criteria is shown in Fig. |||. The distribution of signal events is parameterized with a Crystal Ball 
function (7j, which incorporates a Gaussian core with a power-law tail toward lower masses, and 
accounts for calorimeter shower shape fluctuations and energy leakage. The background is modeled 
by a threshold function 

4.2 Selection of B° -> D^+D*~ Decays 

D* + candidates used in the partial reconstruction of the decay B° — ► D* + D*~ must satisfy |AM — 
AM pcak | < 2.5 (Jam, where AM pcak is the peak of the signal AM distribution observed in the data, 
and a am = 5.7±0.3 MeV/ c 2 is its r.m.s. The CM momentum of the candidate is required to be 

greater than 1.5GeV/c. candidates satisfying these criteria, in addition to those described in 



i 1 (*) I 

section 4.1, are then combined with tt~ candidates to form partially reconstructed B u — > Dl D 



candidates. 

Due to the high combinatorial background in the AM distribution, more than one D* + ir~ 
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candidate pair per event is found in 20% of the events. To select the best candidate in the event, 
the following \ 2 

a _ / Mj-M? e3k \ 2 / M Ds -M%f y / AM — AM peak 

is calculated for each D* + candidate, where M, is the invariant mass of the intermediate 0, K* u , 
or Kg candidate, depending on the Df decay mode, M? eak is the corresponding peak of the signal 
Mi distribution, and Cj is its width. Only the candidate with the smallest value of x 2 in the event 
is accepted. 




5 Results 

The missing mass distributions of partially reconstructed B u — » Ds decays are shown in 

Fig. ||. A clear signal peak is observed in all modes. We perform a binned maximum likelihood 
fit of these distributions. The fit function is the sum of a Gaussian distribution and a background 
function /g given by 

f B (M miss ) = c i(Mo-M miM )^ (3) 
C 3 +(M -M mlss ) 2 

where Cj are parameters determined by the fit, and Mo = M^* — M n = 1.871 GeV/c 2 is the 
kinematic end point. The fits find 3704 ± 232 and 1493 ± 95 peaking events under the Gaussian 
peak in the sum of the D+tt~ and D* + tt~ plots, respectively (Figs. ||,||). However, due to the 
presence of peaking backgrounds, discussed below, further calculation is needed in order to extract 
the signal yields and the branching fractions. 



5.1 Background Study 

We use a Monte Carlo simulation, which includes both BB and qq continuum events, to study the 
missing mass distributions of the backgrounds. We consider two kinds of backgrounds: "peaking" 
background is enhanced under the signal peak at the high end of the missing mass spectrum, and 
"non-peaking" background has a more uniform missing mass distribution. There are two sources 
of the peaking background: 

• Cross Feed (CF): If the soft photon from D* + — > Dfj decay is not reconstructed, B° — > 
D* s + D*~ decays may lead to an enhancement under the signal peak of the D^tt~ missing 
mass spectrum. Similarly, the B° — > DfD*~ decays may lead to a peaking enhancement in 
the D* + ir~ M m ; ss spectrum, due to the combination of a Df with a random photon. 

• Self-Cross Feed (SCF): This is due to true B° -> D* + D*~ decays in which the D+ is 
correctly reconstructed, but combined with a random photon to produce the wrong D* + 
candidate, resulting in a peaking enhancement in the D* + ir~ spectrum. 

Figure || illustrates the missing mass distributions obtained from the Monte Carlo simulation, where 
the cross feed and the self-cross feed are shown separately. Table || presents the reconstruction 
efficiency of correctly reconstructed signal B° — > D*~ decays, as well as cross feed and self- 

cross feed, for events in the signal region M m i ss > 1.86GeV/c 2 . 

In addition to the above backgrounds, we also considered a possible contribution from the 
charged and neutral B decays B — > D s D . These backgrounds were simulated with four D** 
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Figure 3: Missing mass distributions of data events, (a) Dfn with Df — ► 07r + , (b) D* + ir 
with D+ -» </>vr+, (c) D+7T- with £>+ -► (d) D*+tt- with D+ -» (e) D+7T- with 

D+ -» tf^K" 1 -, (f) with L>+ -► The curves show the result of the fit (see text), 

indicating the signal and background contributions. 
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Figure 4: D+ir~ missing mass distributions of 
data events. All the Df decay modes have 
been combined in this plot. The curves show 
the result of the fit (see text) indicating the 
signal and background distributions. 



Figure 5: D* + ir~ missing mass distributions 
of data events. All the decay modes have 
been combined in this plot. The curves show 
the result of the fit (see text) indicating the 
signal and background distributions. 



Table 1: The efficiencies of the partially reconstructed B° — > Di*' + D*~ decays. Columns show the 
contribution of the different generated modes to the D+ir~ and D* + tt~ missing mass distributions 
in the signal region M miss > 1.86GeV/c 2 . Two different B° -> D* s + D*~ Monte Carlo samples have 
been used, one with longitudinal (long.) and the other with transverse (transv.) polarization. 





Reconstructed mode 


Generated mode 


D+7T- 


D*+ir~ 


B u -» DfD*- 


23.6±1.0% 


1.7 ±0.3% 


_» D*+D*- (long.) 


9.0±0.3% 


7.4 ±0.3% 


Self-Cross Feed 




1.6 ±0.1% 


B° -> D* + D*- (transv.) 


10.4±0.3% 


6.9 ±0.3% 


Self-Cross Feed 




1.4 ±0.1% 
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states: D$(j = 1/2), Di(2420), Di(j = 1/2) and L>?j(2460), and their contribution has been 

(*)+ 

determined to be negligible, due mainly to the Dg CM momentum cut. 

Figure |7] shows a comparison of the missing mass distributions in data and Monte Carlo events. 
We assume 1.05% and 1.59% branching fractions for the B° -> D+D*~ and B° -» D*+D*~ decays, 
respectively, in the Monte Carlo simulation. 



5.2 Branching Fractions 

The number of events in the Dfir~ and D* + tt~~ M m j ss peaks is obtained from the fits described in 
section ||. In calculating the branching fractions from these yields, we take into account the fact 
that the peaks consist not only of correctly reconstructed signal, but also of cross feed and self-cross 
feed. This is done by inverting the 2x2 efficiency matrix, whose diagonal elements correspond 
to the sum of signal and self-cross feed efficiencies presented in Table [j], and whose off-diagonal 
terms are the cross-feed efficiencies. The efficiencies corresponding to transverse and longitudinal 
polarization of B° — > D* + D*~ have been weighted according to the measured polarization (see 
section pOj ) . With this procedure, the B° — » D*~ branching fractions are determined to be 

B(B° -► DfD*-) = (1.03 ± 0.14 (stat) ± 0.13 (syst) ± 0.26 (syst B{Df -> (fm + )))% , (4) 

B(B° -» D* + D*-) = (1.97 ± 0.15 (stat) ± 0.30 (syst) ± 0.49 (syst B(D+ -> (fm + )))% (5) 

and their sum is 

UB(B° -» Di* )+ D*~) = (3.00 ± 0.19 (stat) ± 0.39 (syst) ± 0.75 (syst B(Df -> 0tt + )))%, (6) 

where the first error is statistical, the second is the systematic error from all sources other than 
the uncertainty in the — * (pn + branching fraction, and the third error, which is dominant, is 
due the uncertainty in the Df — > <fiir + branching fraction B(Df — > 4>ir + ) = (3.6 ± 0.9)% Q. The 
sources of the systematic error are discussed in section |5.4j . 



5.3 Polarization 

The measurement of the fraction of the longitudinal polarization Tl/T in the B° — > D* + D*~ de- 
cay mode is performed using the events reconstructed in this mode in the signal region (M m i ss > 
1.86GeV/c 2 ). To reduce the systematic error due to large backgrounds, the polarization measure- 
ment is done with only the channel — > (j)Ti + , which has the best signal to background ratio. 
Two angles are used: the helicity angle # 7 between the D*~ and the soft photon direction in the 
D* + rest frame, and the helicity angle 6 n between the D* + and the soft pion direction in the D*~ 
rest frame. Since the B meson is not fully reconstructed, we compute # 7 and 9^ by constraining 
-^miss to the nominal D° mass to obtain a unique kinematical solution for the azimuth (/>. 

The two dimensional distribution (cos0 7 , cos9 n ) is divided in five bins in each dimension. The 
combinatorial background, as well as the cross feed and the self-cross feed obtained using the 
Monte Carlo simulation, are subtracted from this two-dimensional data distribution. The resulting 
signal distribution is corrected bin-by-bin for the detector efficiency, which is obtained from the 
simulation separately for each bin. A two-dimensional binned minimum-^ 2 fit is then performed 
on the efficiency-corrected signal distribution using the fit function 

^ Y± CO s 2 ^sin 2 7 + (l - ^ )sin 2 g 7r 1 + C ° s2 ^ . (7) 

dcos9 n c/cos6» 7 r 7 v r ; 4 w 
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Figure 6: The missing mass distribution of (a) Dfn~ and (b) D* s + n~ Monte Carlo. From top to 
bottom, the overlaid histograms show the contributions of signal, cross feed (CF), self-cross feed 
(SCF) ( only for D* + 7r~) and combinatorial background. 
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Figure 7: The missing mass distribution of (a) Dfir~ and (b) D* + ir~ combinations for data (data 
points) and Monte Carlo (histogram). The contributions from the BB, cc and uds are shown 
separately. The CF and SCF backgrounds are included in the total histogram, not in the hatched 
BB histogram. 
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The resulting fit has a \ 2 °f 23.1 for 25 bins with two floating parameters (Tl/T and total nor- 
malization). Fig || shows the data and the result of the fit projected on the cos# 7 and cosO^ 
axes. 

From the fit, the fraction of a longitudinal polarization is determined to be 

r L /T = (51.9 ±5.0 ±2.8)%, (8) 
where the first error is statistical and the second is systematic. 

5.4 Systematic Errors 



Table 2: Sources of systematic error (%) for B° — > Dg D* branching fractions and B' 
D* + D*~ polarization. 



Source B u - 


-> DfD*~ 


B u -> D* + D*- 




Background subtraction 


2.7 


5.9 


0.5 


Monte Carlo statistics 


4.2 


6.0 


2.7 


Polarization uncertainty 


0.8 


0.5 




Cross Feed 


3.2 


2.4 




% 


1.6 


1.6 




B(4>^K + K-) 


1.6 


1.6 




Particle identification 


1.0 


1.0 


0.1 


Tracking efficiency 


3.6 


3.6 


0.5 


Soft pion efficiency 


2.0 


2.0 


0.2 


Relative branching fractions 


10.2 


10.2 




B(D*+^Dfj) ' 




2.7 




Photon efficiency 




1.3 


0.1 


7T° veto 




2.7 


0.3 


Total systematic error 


13.1 


15.1 


2.8 



The various contributions to the systematic errors of the branching fractions and polarization 
measurement are summarized in Table ^. The dominant systematic error is due to the uncertainty 
in our knowledge of the three Df decay branching fractions. To evaluate the uncertainty due to 
the background subtraction, the signal yield is determined in an alternative way, by counting the 
number of events in the histogram after a bin-by-bin subtraction of the background, determined 
from the Monte Carlo simulation. The difference of the signal yields obtained in this way from the 
results of section |5] is taken as a systematic error. This also accounts for the systematic error due 
to a possible deviation of the signal shape from a Gaussian. 

The Monte Carlo statistical errors in the determination of the signal and the cross feed efficien- 
cies are propagated to the systematic error. The uncertainty in the calculation of the B° — > D* + D*~ 
polarization is propagated to the branching fraction systematic error. The systematic error due to 
charged particle reconstruction efficiency error is 1.2% times the number of charged particles in the 
decay. An additional error of 1.6% is added in quadrature to account for the uncertainty in the 
reconstruction efficiency of the soft pion. 

The systematic error due to the excluding it overlap (7r° veto) requirement was studied by 
measuring the relative yield of inclusive D* + production in data and Monte Carlo events. To 
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Figure 8: Projections of the number of background-subtracted data events on the cos 6? 7 and cos 9. 
axes. The result of the two-dimensional fit is overlaid. 
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evaluate this error, the selection with and without the 7r° veto was applied for the final photon 
from D* s + — > D+j decay. 

For the polarization measurement, the level of the various backgrounds depends on the charged, 
neutral and particle identification efficiencies. The fit was repeated varying the background ac- 
cording to the errors in these efficiencies, and the resulting variations in Tl/T were taken as the 
associated systematic error. 

To check that the simulation accurately reproduces the background M m ; ss distributions in the 
data, a systematic data-Monte Carlo comparison is made in control samples containing no signal 
events. These samples are events with 1.78 < M m ; ss < 1.85GeV/c 2 ; events in the Df sideband 
1.89 < M Ds < 1.95GeV/c 2 or 1.985 < M Ds < 2.05GeV/c 2 ; events in the D* + sideband 170 < 
AM < 300MeV/c 2 ; wrong sign Di*^ + ir + combinations in either the Mp s and AM sidebands or 



signal regions (see section 4.1); and candidates in which M m i ss was calculated using the negative 



(*)+ i 
of the CM D s momentum p The comparison between the data and the Monte Carlo 

simulation of these control samples is shown in Table ||. The discrepancies indicated in Table ||| are 
taken into account in the calculation of the systematic errors. 

Table 3: The average value ((-ZVb — Nmc)/Nmc)j averaged over all bins, where Afo (-/Vmc) is the 
number of data (Monte Carlo) events in a given bin of the M m i ss distribution of the given control 
sample. SB (SR) refers to the Md„ or AM sideband (signal region) control sample. WS indicates 
wrong sign D\ n + combinations, and —p* ( „ )+ indicates that M m ; ss was calculated using the 

negative of the Z?i*^ + CM momentum. The range of the missing mass 1.78 < M m ; ss < 1.87GeV/c 2 
of the control sample is used except for the first line. 

Sample type D+Tr~ D* + ir~ 



1.78 < M miss < 1.85 GeV/c 2 -0.009 ± 0.007 0.075 ± 0.014 

SB -0.075 ± 0.006 0.007 ± 0.022 

SR, WS 0.006 ± 0.008 0.044 ± 0.015 

SB, WS -0.060 ± 0.007 -0.008 ± 0.024 

SR, -p* ( . )+ 0.015 ± 0.009 0.075 ± 0.016 

SB, —p* * M+ -0.062 ±0.007 -0.123 ±0.022 

Average -0.038 ± 0.003 0.032 ± 0.007 



6 Summary 

In summary, using the partial reconstruction technique, we have measured the branching fractions 
B{B° -» D+D*~) = (1.03 ± 0.14 (stat) ± 0.13 (syst) ± 0.26 (syst B{D+ -> </>vr + )))% 

and 

B{B° -» D*+D*-) = (1.97 ± 0.15 (stat) ± 0.30 (syst) ± 0.49 (syst B(Df -» <j)ir + )))%. 
The fraction of the longitudinal D* + polarization in B° — > D* + D*~ is determined to be 

T L /r = (51.9 ± 5.0 (stat) ± 2.8 (syst))%. 
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This measurement is consistent with the theoretical prediction of (53.5±3.3)% § assuming fac- 
torization. Our preliminary results are also in a good agreement with previous experimental re- 
sults |, 0. 
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